Breakdown of universality for unequal-mass Fermi gases with infinite scattering length 
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We treat small trapped unequal-mass two-component Fermi gases at unitarity within a non- 
perturbative microscopic framework and investigate the system properties as functions of the mass 
ratio K, and the numbers A'^i and N2 of heavy and light fermions. While equal-mass Fermi gases 
with infinitely large interspecies s-wave scattering length are universal, we find that unequal 
mass Fermi gases are, for sufficiently large k and in the regime where Efimov physics is absent, 
not universal. In particular, the {Ni,N2) = (2,1) and (3,1) systems exhibit three-body (3b) and 
four-body (4b) resonances a,t k — 12.314(2) and 10.4(2), respectively, as well as surprisingly large 
finite-range (FR) effects. These findings have profound implications for ongoing experimental efforts 
and quantum simulation proposals that utilize unequal-mass atomic Fermi gases. 
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Dilute ultracold atomic Fermi gases have recently at- 
tracted a great deal of attention from the atomic, nu- 
clear, particle and condensed matter communities 
An intriguing aspect of equal-mass two-component Fermi 
gases is their universality. In the regime where the inter- 
species s-wave scattering length Os is much larger than 
the range rp of the underlying two-body (2b) potential, 
the few- and many-body behavior of dilute equal-mass 
Fermi gases is governed by a single microscopic param- 
eter, namely as- Universality has enabled comparisons 
between systems with vastly different length and energy 
scales such as neutron matter, nuclear matter and cold 
atomic gases. Since the interspecies atom-atom s-wave 
scattering length can be tuned to essentially any value 
in the vicinity of a Fano-Feshbach resonance, ultracold 
atomic gases have emerged as an ideal model system with 
which to test fundamental theories such as the BCS-BEC 
crossover theory and to engineer and probe novel quan- 
tum phases [2]. 

A new degree of freedom, the mass ratio k, comes 
into play when one considers unequal-mass Fermi sys- 
tems such as mixtures of non-strange and strange quarks, 
mixtures of electrons with different effective masses 
and dual-species cold atomic gases 0, 0] • For equal num- 
ber of heavy and light fermions, the mismatch of the 
Fermi surfaces gives rise to novel quantum phases such as 
an interior gap superfluid @. Ultracold two-component 
atomic Fermi gases with unequal masses are considered 
to be prominent candidates with which to realize these 
unique phases. Most proposals along these lines assume 
that unequal-mass atomic Fermi systems are stable and 
universal. While this is true for Li-K mixtures, this Let- 
ter shows that these assumptions are, in certain regimes, 
violated. 

We consider few-fermion systems with infinitely large 
interspecies s-wave scattering length a^. The infinitely 
strongly interacting three-fermion system consisting of 
two mass mi fermions and one mass 7712 fermion with 
relative orbital angular momentum L — 1 and parity 
n = — 1 supports, within the zero-range (ZR) frame- 
work, an infinite number of 3b bound states ii k ^ 13.607 



(k — 7711/7712) [3-0I- The properties of these states de- 
pend, as do those of bosonic = O"*" Efimov trimers, on 
Qs and a so-called 3b parameter. For 8.619 < k < 13.607, 
3b resonances have recently been predicted to be acces- 
sible |10i]- We focus on the regime with k < 13.607 
and find: (i) The (2,1) system interacting through a 
purely attractive FR Gaussian interaction potential with 
l/tts = first supports, in the ZR limit, a 3b bound 
state in free space for k = 12.314(2). On resonance, 
the 3b bound state is, similar to a s-wave dimer, in- 
finitely large. Away from the 3b resonance, the behav- 
ior of gas-like states of the trapped system is to a good 
approximation universal, (ii) Adding a light particle to 
the (2, 1) system does not, to within our numerical res- 
olution, lead to a new resonance. The (3, 1) system, in 
contrast, exhibits a 4b resonance at k ~ 10.4(2). On 
the one hand, these few-body resonances open intriguing 
opportunities for studying weakly-bound few-body sys- 
tems. On the other hand, these resonances lead to losses 
in experiments, thereby making the study of macroscopic 
unequal-mass Fermi gases more challenging, (iii) We find 
that unequal-mass few-fermion systems exhibit surpris- 
ingly large FR effects. This finding is relevant since a 
number of numerical techniques are more readily adapted 
to treating FR than ZR interactions. Consequently, a 
full theoretical understanding implies understanding FR 
effects. Furthermore, since realistic atom-atom interac- 
tions have a finite range, comparisons between theory 
and experiment have to account for FR effects. 

Our calculations are performed for a trapped Fermi gas 
with N atoms, N — Ni + N2 {Ni atoms with mass 7771 
and N2 atoms with mass 7772). The model Hamiltonian 
H reads 

^ V 27711 '2 V 

j=l 



j=l k=Ni+l 



where fj denotes the position vector of the jth fermion 
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measured with respect to the trap center and Vtb{rjk) the 



interspecies interaction potential (here, r^-fc 



41). 



Intraspecies interactions are, away from an odd partial 
wave 2b Fano-Feshbach resonance, weak and are ne- 
glected in Eq. ([1]). The spherically symmetric harmonic 
confinement is characterized by the angular frequency 
Lu, which determines the harmonic oscillator length aho, 
aho = \/ft/(2/ia;) with /i = 17111712/ {mi + 1712). Through- 
out, we consider the infinite scattering length limit, i.e., 
l/ttg = 0. Our calculations are performed for the FR po- 
tential Vfr, Vfrir) = ~VQexp[-{r/{V2ro))^], with depth 
Vq (Vq > 0) and range vq {vq ^ Oho)- For fixed rg, we 
adjust Vq so that the 2b system in free space is just at 
the verge of supporting its first s-wave bound state. For 
TV = 3, our results for Vir{r) are compared with those for 
the ZR potential V^^, V^^ir) = 2Tr {fi^ / fi)S{f)-§^r. 

To determine the eigenenergies of the Hamiltonian H, 
we separate off the center-of-mass degrees of freedom 
RcM and solve the Schrodinger equation in the relative 
coordinates. For the unitary system with ZR interac- 
tions, the relative wave function separates into a hy- 
perradial part F^(R) and a hyperangular part ^^(fl), 
= R-^^^~^^/^F^{R)^^{n) [lH; here, R de- 
notes the hyperradius, /ii?^ = J2f=i ^ji^j ~ RcmY-, and 
f2 collectively denotes the remaining 3A^ — 4 degrees of 
freedom. For iV = 3, the eigenvalues of the hyperan- 
gular Schrodinger equation can be obtained by solving 
a transcendental equation (see, e.g., Ref. resulting 

in K(-R) = ^ ^2^1^^^"^^ • = 1~, the quantity Sq — 

defined as the positive root of Sq — decreases from 1.773 to 
as At increases from 1 to 13.607. For k > 13.607, so be- 
comes purely imaginary and Efimov physics emerges 0- 

In a second step, the hyperradial Schrodinger equation 

(if ^ + VAR) + MR) = e^R) (2) 

is solved for F^{R) and E„. This differential equation has 
two linearly independent solutions and gi, |13| . fi/{x) — 
A^x^-+i/2exp(-xV2)ii^i(-g,s^ -t- l,x2) and g^{x) = 
B^x~'^"+^/^ exp(— a;^/2)iFi(— g - s^, —Sy + 1, x^), where 
the non-integer quantum number q is defined through 
El, — {2q + Sy + l)haj and x = i?/(-\/2aho)- We write 
Fu{R) as cos(7r/i^)/iy(i?) — sin(7r/iy)5i,(i?). Requiring that 
F„ (i?) vanishes at large R, the "quantum defect" fj.^, [l3| 
is determined by the condition sin(7r(^y + qj) — 0. Next, 
the allowed values of q are determined by analyzing the 
small R behavior. The function g^, is only normalizable if 
Sy < 1 @,[l3; this implies k > 8.619 for the lowest hyper- 
angular eigenvalue (i.e., v = 0). Thus, for k < 8.619, the 
ZR solution is determined by the exponentially decaying 
piece of fv{R) and the quantization condition gives, in 
agreement with Ref. [ll|, <z = 0, 1, • • • ; we denote the cor- 
responding energy by Ef^i, [solid lines in Fig. (Tfa)]. For 
K > 8.619 (sq < 1), both /j, and g^ can contribute [Tcl| 
and the eigenenergy depends on the boundary condition 
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FIG. 1: (Color online) (2, 1) system with = 1~ at uni- 
tarity. (a) Circles, crosses and pluses show the SV ener- 
gies, extrapolated to the ZR limit, for the three energeti- 
cally lowest-lying states as a function of k. Solid and dot- 
ted lines show Effi [q = Q and 1) and Egfi {q — —so 
and —So + 1), respectively, (b) Pluses, squares, triangles, 
crosses, diamonds and circles show the lowest SV energy for 
K = 10, 12, 12.3, 12.314, 12.4 and 12.5, respectively, as a func- 
tion of ro. Solid lines show fits to the SV energies, (c) Dotted 
and solid lines show l/Lo(a;o) for xq ~ 0.001 as a function of 
q for So = 1/2 (k ^ 12.3131) and So = 11/20 {k ^ 12.0449), 
respectively. The data for so = 11/20 are scaled by a fac- 
tor of 5 for enhanced readibility. In the ZR limit, the reso- 
nance condition implies q — —1/2, 1/2, • • • for so — 1/2 and 
q = -11/20, 9/20, • • • for So = 11/20 (see text for details). 



(BC) at small R. Just as in the case of Efimov trimers, 
this BC is determined by the true atom-atom interactions 
and cannot be derived within the ZR framework. Pa- 
rameterizing the BC by the log derivative Ly{xo), where 
Lu{xq) — [F[,{x) / F„{x)]x=xoj the allowed q values can 
be obtained as a function of L^(xo). Figure [TJc) shows 
[Lo{xo)]~^ as a function of q for sq — 1/2 and 11/20. If 
only gi, contributes, we find, in agreement with note [43] 
of Ref. [11], q — —Si/, —Si, -|- 1, • • • in the ZR limit. The 
corresponding energies Eg^ are shown by dotted lines in 
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Fig. [Ija). In the following, we determine the 3b spec- 
trum for the FR potential Vf^{r) by the stochastic vari- 
ational (SV) method, which makes no assumption about 
the small R behavior of Fi,{R), and interpret it using 
Fig.HTc). 

The SV approach [14] expands the relative wave 
function in terms of a basis set. The proper anti- 
symmetrization of the basis functions, which are opti- 
mized semi-stochastically, is enforced explicitly through 
the application of permutation operators. The resulting 
eigenenergies E{Ni, N2) of the relative Hamiltonian pro- 
vide an upper bound to the exact eigenenergies. The 
functional form of the basis functions used depends on 
the state of interest [3]. For natural parity states, we 
use a spherical harmonic that depends on a generalized 
coordinate and multiply it by a product of Gaussians in 
the relative distance coordinates. In this case, the basis 
functions have definite parity and angular momentum. 
To describe unnatural parity states, we employ so-called 
geminal type basis functions that have neither good par- 
ity nor good angular momentum and select the state of 
interest from the entirety of states. 

Symbols in Fig. [TJb) show selected SV energies for the 
(2,1) system with = 1~ at unitarity as a function 
of ro. To extrapolate to the ro — limit, we per- 
form 4-5 parameter fits to the SV energies [solid lines 
in Fig. mb)]. The resulting extrapolated ZR energies 
are shown by symbols in Fig. [T][a). The key characteris- 
tics of Figs. [TJa) andlljb) can be summarized as follows: 
(i) The dependence of the 3b energy on ro increases as 
K increases from 1 to about 12. For k = 12, e.g., the 
difference between the extrapolated ZR energy and the 
energy for ro = O.Oloho is about 30%. For k = 12.314, 
the energy depends comparatively weakly on ro . (ii) For 
K — 12.4 and 12.5, the ground state energy is negative for 
sufficiently small ro and diverges as r^'^ with decreasing 
ro. (iii) For k 7^ 12.314 and k — 12.314, the extrapo- 
lated ZR energies [symbols in Fig. [Ija)] agree to a very 
good approximation with Ef^ and Eg^o, respectively. We 
interpret the dropping of the energy around k « 12.3 
(so ~ 1/2) as a 3b resonance. The resonance position is 
found to be K = 12.314(2). 

Although the occurance of the 3b resonance depends 
on the underlying 2b potential, we now argue that a 3b 
resonance occurs most likely if sq ~ 1/2. The hyperradial 
wave function of the 3b system with FR interactions is for 
i? > ro given by F^. Figure[TJc) shows the dependence of 
the allowed quantum numbers q on the value of l/Lo(a;o) 
for xq = 0.001. For sq ^ 1/2 [the solid line in Fig. [Tic) 
shows an example for sq = 11/20], nearly all values of 
the logderivative result in q ss 0, 1, • • • , implying that 
the realization of a 3b resonance for sq ^ 1/2 requires 
careful fine-tuning (as xq decreases, the "corners" of the 
solid line near integer q values become even sharper). If 
realized, such a resonance is narrow, with the size of the 
zero-energy trimer set — as in the case of non-universal 
p-wave dimers — by the effective angular momentum bar- 
rier. Physically, an unnaturally large contribution of gi, 
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FIG. 2: (Color online) Squares, triangles and diamonds show 
the extrapolated ZR ground state energies of the (2, 2) system 
with = 0^, 1~ and 2"*", respectively. The extrapolation to 
the ZR limit is based on 3-5 parameter fits to the SV energies 
for ro ~ O.OOSflho to 0.02aho. Circles and a soUd line show the 
energies of the (2, 1) system with = 1~ [see Fig. [Ha)]. 



is, for So > 1/2, "suppressed" by the effective repulsive 
angular momentum barrier in the hyperradial coordinate. 
For So — 1/2 [dotted line in Fig. [TJc)], in contrast, the 
resonance is broad and on resonance the size of the zero- 
energy trimer in free space is, much like that of s-wave 
dimers at unitarity, infinite. Our analysis of the struc- 
tural properties such as the hyperradial density for FR 
systems confirms this conclusion. Since the effective an- 
gular momentum barrier vanishes for sq = 1/2, is 
no longer naturally suppressed and the existence of a 3b 
resonane is more probable. While our analysis shows 
a 3b resonance occurs more likely if so ~ 1/2, the ZR 
model cannot predict whether or not such a resonance 
does indeed occur. Based on the FR results presented 
in Figs. [Ija) and (b), we speculate that other classes of 
interaction potentials likely also support a 3b resonance 
near sq ~ 1/2. 

We now investigate the (2,2) and (3,1) systems, i.e., 
we add respectively a light and a heavy atom to the (2, 1) 
system. Symbols in Fig. [H^a) show the extrapolated ZR 
energies for the lowest natural parity states of the (2, 2) 
system with L = — 2. The FR effects are comparable 
to those of the (2,1) system discussed in the context of 
Fig. [1] . The uncertainty of the S V energies is less than 
1% and the uncertainty of the extrapolated ZR energies 
is primarily due to the fact that our SV calculations for 
the (2,2) systems are limited to ro > 0.005aho- For all 
three L considered, the (2, 2) energies lie above the (2, 1) 
energies and show a notable drop at k « 12.3, which 
we attribute to the presence of the 3b resonance. Our 
calculations show no evidence for a (2,2) resonance. We 
also considered unnatural states of the (2, 2) system and 
find that the lowest unnatural parity state lies above the 
lowest natural parity state. 

We now show that the (3, 1) system exhibits a 4b res- 
onance at a K value that differs from that at which the 
3b resonance occurs. The energetically lowest lying state 
of the (3, 1) system has L = 1 and unnatural parity. Our 
extrapolated ZR energies are 5.1huj, A.Qhu, 4.5haj and 
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3.9huj for K = 1, 2, 4 and 8, respectively. While our 
calculations for larger k do not allow for a reliable ex- 
trapolation to the ZR limit, they do allow for a reliable 
determination of the (3, 1) resonance position. To locate 
the resonance position, we monitor whether the SV ener- 
gies decrease or increase with decreasing tq. For k = 10, 
e.g., we find that the energy increases with decreasing 
tq. For K = 10.6, in contrast, we find negative energies 
for ^0 ^5 0.02aho. Performing additional calculations for 
K = 10.3, 10.4 and 10.5, we find that the (3, 1) resonance 
is located at k = 10.4(2). As in the (2, 1) case, the exact 
position of the (3, 1) resonance depends on the details of 
the underlying 2b interactions. In particular, the analy- 
sis of the small R BC of the hyperradial wave function of 
the (3,1) system parallels that of the (2,1) system. This 
implies that the smallest mass ratio at which a (3,1) 
resonance can occur depends on the solution to the hy- 
perangular Schrodinger equation, i.e., the sq of the (3, 1) 
system needs to be smaller than 1. It appears likely that 
the (4, 1), (5, 1), etc. systems exhibit resonances at suc- 
cessively smaller k. 

In summary, we considered unequal-mass two- 
component Fermi gases interacting through FR poten- 
tials with infinite interspecies s-wave scattering length 
and found 3b and 4b resonances in regimes where Efimov 
physics is absent. These resonances are non-universal in 
the sense that their exact position and properties depend 



on, besides a^, at least one additional parameter. We ar- 
gued that 3b resonances occur most likely if k « 12.3 
(sq ~ 1/2) and that 4b systems likely exhibit, assum- 
ing a 3b resonance exists near sq « 1/2, a 4b resonance 
near k « 10.4. While we adjusted sq by varying k, ex- 
perimentalists could tune the A^-body sq by utilizing an 
intraspecies Fano-Feshbach resonance jT(l |. Importantly, 
our results for the (A^i,l) systems apply not only to 
Fermi-Fermi mixtures but also to Fermi-Bose mixtures, 
making the ^Li-^^Sr system (k « 12.4) — thanks to the ex- 
istence of optical Sr-Sr Fano-Feshbach resonances [l^ — a 
promising candidate for experimental studies. We note 
that the 3b resonance discussed here differs from the res- 
onances discussed in Ref. [l^ for positive Og. The break- 
down of universality and the instability of unequal-mass 
atomic Fermi gases discussed in this work also provides 
opportunities: Our theoretical study paves the way for 
exciting investigations of novel few- and many-body sys- 
tems with simultaneous 2b and 3b or 2b and 4b reso- 
nances, either in a dipole trap or an optical lattice set-up. 

After completion of our manuscript, we received a 
related, independent manuscript by S. Gandolfi and J. 
Carlson [l^; see also talk by J. Carlson on 03/11/10, 
http:/ / www. int. Washington, edu/ talks /Workshops / int_ 1 0-4 6 W/. 
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